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We derive a differential equation for the one-particle-irreducible vertex functions of interacting 
fermions as a function of the temperature. Formally, these equations correspond to a Wilsonian 
renormalization group (RG) scheme which uses the temperature as an explicit scale parameter. Our 
novel method allows us to analyze the competition between superconducting and various magnetic 
Fermi surface instabilities in the one-loop approximation. In particular this includes ferromagnetic 
fluctuations, which are difficult to treat on an equal footing in conventional Wilsonian momentum 
space techniques. Applying the scheme to the two-dimensional t-t' Hubbard model we investigate 
the RG flow of the interactions at the van Hove filling with varying next-nearest neighbor hopping t' . 
Starting at t' = we describe the evolution of the flow to strong coupling from an antiferromagnetic 
nesting regime over a d-wave regime at moderate t' to a ferromagnetic region at larger absolute 
values of t' . Upon increasing the particle density in the latter regime the ferromagnetic tendencies 
are cut off and the leading instability occurs in the triplet superconducting pairing channel. 



I. INTRODUCTION 

Wilsonian momentum shell renormalization group (RG) schemes have become an important tool for the weak- 
coupling analysis of interacting electron systems because they provide a means to sum competing classes of diagrams 
in perturbation theory in a consistent manner. Moreover, with the Wilsonian RG schemes one can systematically 
approach potential singularities, which would cause problems in a direct perturbation expansion. Successes of RG 
methods in interacting electron systems are for example clear arguments for the stability of the Landau-Fermi liquid 
in two and more spatial dimensions for short range repulsivei-jnteractions and non-ncstcd Fermi surfaces unless 
a Kohn-Luttinger instability intervenes at lowest temperatureaa 13, or the classification of ground states in quasi- 
one-dimensional systemgcl U. Moreover, numerous one-loop RG calculations have been used recently to analyze the 
interplay between antifierromagnetic (AF) and d-wave superconducting (SC) pairing tendencies in the two-dimensional 
(2D) Hubbard modelQ tj. On the other hand, the Hubbard model was originally introduced to describe ferromagnetism 
(FM) of itinerant ekctrons (see, e.g., Refs. p^,p^). The existence of spin-polarized ground states has been proven 



in one-dimensioiiajti'E3 and in special higher-dimensionalEj models and also established in the limit of infinitely 
many dipiension£ln3. In the t — t' Hubbard model on the 2D square lattice at weak-ta moderate U, Hartree-Fock 
resultaHj, T-matrix approximation together with Quantum Montecarlo sirpulationsoEj and a generalized random 
phase approximation approach including local particle-particle correlationsc3 point towards a FM regime for larger 
absolute values of t' around the van Hove-BJling where the Fermi surface (FS) contains the saddle points at (±7r,0) 
and (0, ±7r). Recent parquet calculations^ indicated similar tendencies. However, momentum-shell RG methods 
have up to now not produced any evidence confirming these findings. 

In this paper we introduce a novel RG scheme, which we call the temperature-flow RG. Its distinguishing feature 
. ^ is that the temperature itself plays the role of the flow parameter. Indeed, it has been observed in many studies that 
. varying a cutoff scale is in many ways similar to varying the temperature, but here we provide an exact RG flow 
" equation for the one-particle-irreducible correlation functions, parameterized by the temperature. Using this method, 
, we show that in the one-loop approximation, one finds a ferromagnetic regime in the t — t' Hubbard model. More 
precisely, the new flow is similar to previous ones for small t' and near to half-flUing, but for larger \t'\, a ferromagnetic 
regime appears. The boundary of the FM phase is roughly at the same t'^ where it would be in the Hartree-Fock 
approximation, but with our RG calculation we find a quantum critical point separating the FM regime from a d-wave 
superconducting regime instead of a first order transition to an AF phase as in Hartree-Fock. We also find a triplet 
superconducting state when the filling is increased above the van Hove value. 

In the following, we discuss the reasons why this FM phase was not seen in previous RG calculations. The main 
point is that when the flow parameter is given by an infrared cutoff, the particle-hole excitations with small total 
momentum are artificially suppressed above the scale set by the temperature, but in many interesting cases one cannot 
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run the flow down to these low scales because the one-loop RG flow of the interactions leads to strong coupling and 
even diverges at a nonzero scale. 

In RG flows to strong coupling, the main interest is to determine the nature of the strongly coupled state. Very often 
there are several candidates for the latter and it is natural to assume that the strong coupling state with the highest 
energy gain will be realized. In the RG treatment this state is typically associated with the dominant channel at the 
energy scale at which the RG flow diverges, e.g. determined through an analysis of the corresponding susceptibilities. 
However, one-loop flows become unreliable once the coupling constants get larger than a value which depends, among 
other factors, on the Fermi surface. Thus it can happen that the one-loop approximation breaks down before processes 
involving particle-hole fluctuations at small wave-vectors, which could drive ferromagnetism or screen longer-ranged 
interactions, can intervene. Because of its inequivalent treatment of FM fluctuations on the one side and AF and SC 
tendencies on the other side, this type of momentum shell RG with infrared cutoff does presumably not represent, |-aJb 
least when used too straightforwardly, an appropriate method to answer the question which of the channels prevailsE^. 

Thus, although momentum shell schemes apply, and can even be controlled rigorously, as long as the running 
coupling constants do not get large, it is interesting to investigate and compare alternative schemes which are not 
based on strict infrared cutoffs and therefore capable of taking into account low-energy fluctuations at all wave- vectors 
in the same way. 

We now speciflcally discuss the Stoner problem. Consider an electron system with a density of states p per spin 
orientation and a repulsive onsite interaction [/ > 0. The RPA spin susceptibility is 



where Xs(a;, q) is the bare particle-hole bubble 



d'^k nF{e{k + q}) - nF{e{k)) 
(2^ lu; + e{k + q) - e{k) 



with np the Fermi distribution. Sending w — s- and afterwards q—>0, we get 

X°(a> = 0,<f^O) = - / ^lnF{e{k))= f dEp[E)5T{E) (3) 



(27r)^ dt 

where St{E) — -~^np{E) — (4Tcosh^(i?/2r))^^ becomes the delta function as T ^ 0. If the density of states p is 
continuous, we get, at T = 0, XsC"^ = 0, (f ^ 0) = p{Q)/{l — p{0)U). The pole in this expression leads to the simple 
Stoner criterion 

Up{0) > 1 (4) 

for a ferromagnetic ground state. 

Of course, it is well known that the Stoner criterion overestimates the tendency towards ferromagnetism, and one 
should include correlation effects. One of them is that in (^, U is expected to get replaced by an_ effective Stoner 
coupling C/stoncr- This was first discussed in the so-called T-matrix approximation by KanamoriEJ. On the other 
hand, the T-matrix approximation also sums only a subclass of diagrams of one-loop type. Thus the natural question 
arises whether by RG methods one could perform a consistent resummation of all one-loop contributions, to treat 
the ferromagnetic tendencies and the screening together with other possibly competing tendencies. 

Let us attempt to apply a one-loop momentum-shell scheme to the Stoner problem. The divergence of the uni- 
form spin susceptibility is equivalent to a pole in the RPA perturbation series for the forward scattering amplitude 
a°(fc, k' , q-^0) which is obtained from the antisymmetric irreducible two-particle scattering vertex r*-^-* via 

a"(fc, k') = lim r('') \{k + q, s, k' , ~s) (fc, -s, k' + q, s)] (5) 

= lim (r('^) \{k + q, s, k', s) {k, s, k' + q, s)] - F^''^ \{k + q, s, fc', -s) {k, s, k' + q, -s)j \ . 
ij-^o L L J L J J 

Since the (f — > particle- hole loop is the only diagram occuring in the RPA calculation of a°-{k, k') and also in the 
RPA calculation leading to the Stoner criterion, we can limit the RG analysis to this channel, in complete analogy to 
Cooper or SDW (spin density wave) instabilities, where the ladder summation of the zero-total-momentum particle- 
particle diagram or the particle-hole diagram with the appropriate momentum transfer, respectively, are written as 
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solutions of the one-loop RG equation for the corresponding channel only (as indicated above, we do a study taking 
into account all contributions below). The RG equation is then the differential equation 

I^C/a = -R^Ul (6) 

where R\(uj,q) is given by the integral for Xs? but with the integrand containing an additional factor ■^[h{e{k + 
q) / A)h{e{k) / A)], with h{x) the cutoff function, say, h{x) ~ for |a;| < 1. In the limit of interest for FM, 

Ra{u; = 0,q^ 0) - ^^A(r) , i^A(T) = J dEp{E)5T{E) h{E/Af , (7) 

and (^) integrates to 

= l-C/A„(FA(r)-FA„(T)) 

If Aq is larger than the bandwidth, F/^g{T) = 0. For A ^ 0, F\{T) Xs(0, 0), given by (|3|), so we recover exactly 
the RPA result, and thus the Stoner criterion. 

However, a closer look at reveals that Fa{T) is negligibly small for all scales A above 0{T) because 6t{E) < 
^-\E\/T g^^^ ^YiQ cutoff function h restricts to \E\ > A. Thus the contribution to the RPA is built up only at the very 
lowest scales pfrilie flow: small-g particlc-holc fluctuations, such as ferromagnetic tendencies, occur only at low RG 
scales A w TOt^. For the derivation of the RPA as a solution to an approximate RG, this is not a problem. But 
our goal is to understand the interplay of different instabilities, and it can happen that another instability driven by 
terms that contribute to the flow at higher scales gives rise to a flow to strong coupling at scales A > T. Then the 
small-q contributions would remain unintegrated because, as discussed above, we have to stop the one-loop scheme 
once the couplings cease to be small. This appears to be particularly unjustified if the density of states (DOS) at the 
FS is much larger than the DOS at higher band energies, e.g. in cases where the FS is close to a van Hove singularity. 

One can argue that this problem is not only due to the one-loop approximation but that it also has to do with 
restricting to g = 0. If \q\ is of order A/vp, the infrared cutoff does not suppress the particle-hole susceptibility any 
more. Thus one problem is the nonuniformity in the buildup of the function of q, and one might think that averaging 
over 5 in a region of size 0{A/vf) may help. 

However, in this paper we attack the problem directly: as we shall show, in the temperature-flow RG scheme, 
the small q fluctuations are properly taken into account directly, without further approximations like the averaging 
mentioned above. 

In the next section, we set up the exact RG equation for the temperature-flow, and describe the one-loop approxi- 
mation. After considering simple examples we apply the scheme to the two-dimensional t-t' repulsive Hubbard model 
on the square lattice. This model is interesting because it supports a variety of fluctuations which can become singular 
at low temperatures and an unbiased treatment is highly desirable. On one hand it shows strong antiferromagnetic 
fluctuations due to partial nesting of the Fermi surface at finite wavevectors and omnipresent superconducting insta- 
bilities at low temperatures are enhanced by the anisotropic scattering. On the other hand the 2D Hubbard model is 
also a good candidate for a ferromagnetic ground state even at weak coupling because of the large density of states 
at the Fermi surface for certain particle densities. 

II. FORMALISM 

Here we derive differential equations which determine the evolution of the IPI (one-particle-irreducible) vertex 
functions with varying temperature. A perturbative RG scheme describing the change of the IPI functions with 
variation of a general parameter-dependent quadratic part of the free action, s'^^ , was described in Refs. H] and [l^ . 

(2) 

In these applications the continuous change in corresponded to integrating out modes with band energy at the 
running RG scale A. Here we show how the same RG equations can be used to describe the variation of the IPI vertex 
functions with changing temperature. In this case the scale parameter can be directly taken as the temperature of 
the model system. 

Consider a fermionic system on a 2D lattice at temperature T, with action 

-/2 



/d k - ^ ^ 
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ki, s — ^ 



k2,s' 



FIG. 1. The coupling function V{ki,k2,k3). ki is fixed by wavevector conservation on the lattice. 



■4 , «w„4 )ips'{k2, iw„2 ) V's (fci , iw„ J 



J2 



/d k [ - ^ ^ ^ 



(9) 



Here ■0 and are anticommuting Grassmann fields and the Matsubara frequencies a;„ = (2n + l)7rr are summed 
over all integer n. We have assumed spin-rotation-invariant and frequency-independent interactions. Further we 
have added source fields ^ and ^. Obviously both the quadratic and the quartic part of the action depend on the 
temperature T. 

We now recall the case where only the quadratic part of the action depends on an additional continuous parameter 
A. Formally we go from S" to a A-dependent action by replacing 



T 



iLO — e{k) 



(3a(zcj, k). 



(10) 



One can then derive RG differential equations describing the variation of the IPI vertex functions with A. These 
equations are obtained from the A-dependent Legendre transform 



rA(</.>) = TyA(^,0-EE 



<pk 

(2^ 



(t>s{k,iujn)^s{k,iujn) -f (f>s(k,iujn)is{k,iujn) 



of the generating functional of the connected non-amputated m-point Green functions W'^"^\ 



(11) 



(12) 



Here we abbreviate the notation and define a combined index p — {iojn, k, s), moreover J dp = ^„ ^ / -^fjryr- The 
RG differential equations describing the evolution of the selfenergy Sa(p) and coupling function Va(pi,P2,P3) with 
Pi = Pi + P2 ^ P3 (from which the fully antisymmetric four-point vertex can be reconstructed^, see also Fig. ||) with 
variation of that parameter read 



d f 

— Ea(p)= j dp'SAip') [VA{p,p',p')-2VAip,p\p)] 



and 



with 



dA 



VAipi,P2,P3) = Tpp^A + r^HA + '^PH,. 



(13) 



(14) 



TpPA{pi,P2]P3,Pi) = 

dpVAipi,P2,p) L{p, -p + pi +P2) Va{p, -p + Pl +P2,P3) 



(15) 
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FIG. 2. The particle-particle and particle-hole diagrams contributing to the one-loop RG equation. 



'^PH,a{Pi,P2;P3,P4) 



dp 



-2Va{pi,p,P3) L{p,p + pi -P3)Va{p + pi -p3,P2,p) 



+Va{pi,P,P + Pi -p3)L{p,p + pi ~P3)Va{p + pi -p3,p2,p) 
+Va{pi,P,P3)L{p,P + Pi -P3)Va{p2,P + Pi - P3,p) 

^ff,A(Pl'-P2;P3,P4) = 

- j dpVA{pi,P + P2 -p3,p)L{p,p + p2 ~P3)Va{p,P2,P3) 

In these equations, 

L{p,p') = Sa{p)w'^a\p') + W^^\p)Sa{p') 
with the so-called single-scale propagator 



Sa{p) = -W'}^\p) 



dK 



.ip) 



W^'\p). 



(16) 



(17) 



(18) 



(19) 



The one- loop diagrams corresponding to the terms (|l5|), ([l6| ) and (17) are shown in Fig. ^ 

In the typical momentum-shell RG, the varying parameter A is an energy scale which separates high and low 
energy modes. The strategy is to integrate out the high energy modes first. In this case the scale parameter only 
affects the quadratic part Qa{p)i which is multiplied with an appropriate cutoff function. Here, for the reasons 
discussed in the introduction, we want to treat the temperature as varying parameter. Our reasoning is as follows: 
at high temperatures, where ttT is larger than the bandwidth and the interaction energies, perturbation theory 
converges, and moreover the corrections to the selfenergy and four-point function are of order 1/T, hence small. 
Thus at high temperature, the vertex functions are essentially identical to the terms in the action. Then we track 
the renormalization of the vertex functions when the temperature is lowered. Of course this idea is implicit in most 
of the well-known scaling approachesu. On a technical level however this strategy has usually been cast into some 
cutoff- variation procedure with similar results as the more elaborate modern Wilsonian schemes. 

In order to apply the RG scheme for the IPI vertex functions with T as a flow parameter we first have to perform 
a transformation which shifts all temperature dependence to the quadratic part of the action. 



A. New fields 



by 



The T^-factor in the interaction part can be removed by transforming the action to the new fermionic fields given 
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The source fields are chosen to transform according to 
The action then reads 

(27r)2 
d^ki <Pk2 d'^ks 



(20) 
(21) 



n s ^ ' 

■EE 



ni, 712,713 S,s' 



(27r)2 



(27r)2 (27r)2 (27r)2 
fjs{k,n)(s{k,n) + ris{k,n)C,s{k,n) 



V{ki,k2, kz)r]s{kz,nz)ris'{ki,ni)r]s'{k2,n2)-qs{ki,ni) 



(22) 



Only the quadratic part depends on the temperature. The connected m-point correlation functions wjj^"^ and W^^^ 
are related as follows (the change in the integration measure drops out): 

) = {lpsi{kl,iuJnJ ■ . . V's„(fcm,«^^n™)) 



T-3™/'*(?7,i(fci,ni)...?7,„(fc™,?^™)) 
T-'"'/^Wj,"'\k,,n,,...km,nm). 



(23) 



Using the action in terms of the new fields we can now apply the RG formalism from Ref. ^ This yields the flow of 
vertex functions T^\ki,ni, . . .km,nm) when we change the temperature. Since they are the expansion coefficients 
of r^, which is a functional of of fields conjugated to the C source fields, the F^™-* are related to the vertex functions 
of the ip fields via 



in particular 



and 



^ (fci , iujni J ■ ■ ■ k„i , iuJum ) — ^ ^ (fci , rii , . . . kjn , rim) ■, 



ri {ki , iuJni , k2 , iuJn2 ) = T^/^r*^) (fci , rii , 4 , «2 ) 



ri*^(fci, zwnj, . . . ^4, iw„^) = r^r^'*)(fci, ni, . . . fc4, ^4) , 



Let us consider an example. If the interaction term at temperature T is taken to be 



d^ki d'^k2 d^k^ 
(2^ (27r)2 {2tt)' 



■ Vt (fcl , ^2 , fca) V's (fcs , ii^na )'>Ps' (^4 , iuJni)'4's' (fe , «t^«2 ) V's (^1 , i^m ) , 



the antisymmetric four-point vertex in the i/i-fields readsH, to first order in the interaction, 



.(4) 



(/ci, s,zw„i;/c2, s',?a;„2; fc3,s,ia;„3; fc4,s',ia;„J = VT(fci, ^2,^3) - <5s,s'VT(fc2, ^1,^3) 



The four-point vertex for the 77-fields is simply 

r|,'^^(fci, s, ni;k2,s', n2\ h, s, n^; k^, s' , n^) = r^^^j,{ki,s, zw„j; fc2, s', «w„2; ^3, s, iujns^h, s' , iw„J , 



(24) 
(25) 
(26) 

(27) 

(28) 
(29) 



hence of order T*', and it contains the coupling function Vt- Thus the RG equations for the four-point-vertex in terms 
of the 7y-fields exactly yield the flow of the coupling function, and we can use the relations ( p^ to express the results 
in terms of the original vertex functions. 

In this paper, like in Refs. [| and |l|, we will neglect selfenergy corrections to the flow and furthermore truncate 
the infinite hierarchy of RG equations for the m-point vertex functions by setting all IPI vertex function higher than 
fourth order equal to zero. These approximations restrict the validity of the results to the weakly coupled and weakly 
correlated regime. We also drop the frequency dependence of the four-point vertex. 
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B. The single-scale propagator 



The quadratic part of the action expressed in the 7/-fields reads Qn,Tik,n) — T ^/^[iuin — e^:)- In absence of 
selfenergy corrections, we obtain 



(30) 



Therefore the single-scale propagator is given by 



STik,n) = ~W^^lik,n) 



dT 



Qn.xik.n) 



W!^^{k,n) 



iw„ + e(fc) 



iiOn - e(fc) 



(31) 



C. One-loop diagrams 

Next, consider the product of two fermion propagators L{p,p') defined in (p^, which appears in the RG equation 
described in The index p contains the frequency lu and wavevcctor k, at band energy e associated to one propagator 
and p' the frequency oj', wavevector k' at band energy e' associated to the other one. Inserting ( ^0| ) and (|3^) into 
(ITsI), we have 



L{p,p) 



iuj + e 



1 



1 



{iuj — e)^ iuj' — e' iuj — e iiuj' — e'Y 
The terms on the right hand side of the flow equation are of the form 

1 



dT 



T- 



1 



1 



dp L{p, p')<I>(fc, s, fc', s') 



d^k \ - d 

(27r)2 ^dT 



iuJn ~ e{k) iujn' - e(fc') 



$(fc,5,fc',s'), 



(32) 



(33) 



i.e. they are temperature derivatives of bubble diagrams. Thus, in order to calculate the temperature flow of the 
coupling function, we can think in terms of the usual one-loop diagrams of the perturbation expansion in the original 
V'-fields and simply take the temperature derivative of the particle-hole and particle-particle bubbles. 

III. SIMPLE EXAMPLES 

Here we apply the scheme to two simple examples and show that it yields the results expected from direct per- 
turbation expansions. In particular we describe how the temperature-flow scheme allows us to describe Cooper and 
ferromagnetic instabilities in an analogous fashion. 



A. Cooper instability 



Here we show that for a small attractive electron-electron interaction at temperature Tq the formalism leads to a 
Cooper instability at lower temperature Tc- We assume that at T = To the action is given by the BCS model. 



^To^ J tp{k, iuJn) (iuJn 



VF{k - kF))'4'{k,iU!n) 



(34) 



Here, ]/;^°°p'"' < is the attractive interaction between electrons with opposite wavevectors and the bandwidth W is 
taken such that vpkp — W > Q. We assume a linear dispersion. We are interested in the low-temperature behavior 
of the four-point vertex with zero incoming total momentum and frequency, 
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^(4) 
- tp.T 



(k, s, -k, -s) ^ {k', s, -k', -s) = T^Vrik, -k, k') 



Focusing on the pairing channel, we neglect all particle-hole contributions. Therefore we are left with the temperature 
derivative of the particle-particle bubble and Vrik, —k,k') — 1/"^°°^°" remains independent of k and k' . The RG 
equation for Vrp ° °p"' involves the particle-particle bubble n(a;, with lo — and pair momentum q^O. Using the 
results from Sec. [I C it reads 



Cooper 



^Cooper 2 _^n(0, 0) 



We can immediately integrate it, similarly to (0), to get 

1 



Cooper 2 



d 



1 ~ 2nF{e) 
Ye 



Cooper 



To 



W 



2e 



(35) 



(36) 



Assuming that the bandwidth W is much larger than the temperature and that the density of states is a constant poi 
the integral is logarithmic and 



Cooper 



Cooper 



i-Po^^rP"iog(r/ro) 



(37) 



Thus, starting with an attractive interaction V^°°^" < at a higher temperature To ^ we obtain a pole in the 
effective s-wave Cooper pair scattering at 



Tc = To exp 



1 



Po 



V 



Cooper 



To 



(38) 



This result is analogous to the ladder summation from straightforward perturbation theory when we replace Tq 
with the Debye frequency and identify the initial interaction at this temperature with the bare attractive 

electron-electron interaction y^°°p°' ^ e.g. mediated through phonon-exchange. Although these identifications cannot 
be justified directly within a microscopic model with a specific Hamiltonian, they appear to be highly plausible for 
an effective description of the system. 
In fact one can obtain the BCS relation 



Tc = l.IWexp 



1 



It rCoopcr I 

Po\Vto I 



(39) 



as well if one assumes that Tq ^ W. In that case, inserting Tq in the right h;uid side of ( |36| ) gives zero because at 
high temperatures np ~ 1/2. Thus we recover the linearized BCS gap equationEj at T^,, 



-w 



2e 



(40) 



as a condition for a divergence of y^°°P°' at T = Tc. From solving (EH), one obtains M). 



B. Stoner instability 

Next we consider the instability towards ferromagnetism in a Hubbard-type system with repulsive onsite interaction 
UcS- We show that for a logarithmically divergent density of states within an effective bandwidth WcS, 

p(6) = p-log^, (41) 

an approximation to the RG equation reproduces the Stoner instability obtained by straightforward perturbation 
theory. In the latter treatment one sums up the particle-hole bubbles for momentum transfer q ^ with the 
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intermediate particles having band energies within ±Wcff. The van Hove singularity in the density of states causes a 
breakdown of the paramagnetic state at temperatures below 



Tc ~ WcS exp 



UcSP 



(42) 



Analogous to the diagram summation we restrict the RG analysis to the crossed particle-hole channel with small 
momentum transfer. The Stoner instability should then appear as a divergence of the forward scattering amplitude 
a°(A:, fc', (f— > 0), which is in our RG language given by 



lim r^'^^ 

q— »0 ^ ' 



ik, s, k' + q, -s) (k', s,k + q, -s) = lim T^Vrik, k' + q, k') . 



Next we will assume an isotropic system and write 



Yiu\VT{k,k' + q,k') 



V™. 



Using Eq. (^) with the RG equation for this coupling function reads 



dT ^ 



d 

/ de p{e) —npie) . 
J-W de 



(43) 



Note the similarity of, but also the difference between, ( ^3|) and (^ - both equations have the same form, but there 
is no cutoff function in (B^) . Integration of (|43|) gives 



= F{T) - F{To), 



with 



Because 6t{E) oc T ^, F{To) vanishes as Tq oo, hence we recover for large Tq the RPA result 



F{T)^X°{0,0)= / dEp{E)5T{E) 



1 - V™F{T) 



(44) 



(45) 



(46) 



In the case of the density of states (^), we can also use integration by parts to get, for Tq -C W, 



d 

dT 



T^FM 2- 



J-w ^ dl 1 J_w de 



The resulting differential equation is solved by 



T/FM 



1 



y™piog(r/ro) 



(47) 



Analogous to the Cooper instability studied above, the initial interaction Vtq at higher temperature Tq <C , in this 
case repulsive, leads to a pole at 



Tc = Tq exp 



1 



pV^^ 



'To 



(48) 



signaling the instability of the paramagnetic state with respect to spontaneous polarization. Again this result is in 
close analogy with the RPA summation when we replace the initial temperature Tq with the effective bandwidth 
Weff. Of course this approximation to the RG equation can also lead to pole for a non-divergent density of states if 
the initial interaction is sufficiently strong. In this case, upon integrating from Tq W down towards T = 0, we 
encounter a pole if the simple Stoner criterion (^ is fulfilled. 

Thus, in our RG treatment, the simple Stoner criterion can be derived as a consequence of an approximation to 
the one-loop RG equation where only the crossed particle-hole terms are kept. As we shall see in the next section, 
the full one-loop flow is more complicated, and one can thereby see the limitations of the naive Stoner criterion, and 
instances where it becomes invalid because of superconducting and antiferromagnetic tendencies. 
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FIG. 3. The Brillouin zone, Fermi (solid line) and umklapp (dashed line) surface and the lines in the patch lines for TV = 32. 
The circles denote the wave vectors for which the four-point vertex function is calculated. For all wave vectors inside the same 
patch (shaded region) around the lines, VA{ki, k2, ka) is approximated by a constant. 

IV. APPLICATION TO THE 2D HUBBARD MODEL 

Here we apply the RG scheme to the 2D t-t' Hubbard model. In our setup with the temperature as a flow parameter, 
this again means that we assume that at some higher temperature Tq the single-particle Green's function of the system 
is adequately described by 

Go(zo;,fc) = ^—r- (49) 

iuj — e(fc) 

with the band dispersion e(fc) = —2t [cos + cos ky] — At' cos k^ cos ky — ^ including nearest and next-nearest neighbor 
hoppings t and t' , respectively. The scattering vertex is given by a local repulsion 

VT„iki,k2,k3) . (50) 

Then by following the flow of the scattering vertex with lowering the temperature, we obtain information on the 
possible low-temperature phases of the system. In particular we analyze which classes of coupling functions and 
which susceptibilities become important at low temperatures. For a large parameter range, we observe a flow to strong 
coupling. This means that at sufficiently low temperatures some components of the coupling function Vriki, ^2, k^) 
reach values larger than the bandwidth « 8t. The approximations made, i.e. the restriction to one-loop equations and 
the neglect of selfenergy corrections, fail when the couplings get too large. Therefore we stop the flow when the largest 
coupling exceeds a high value larger than the bandwidth, e.g. Vr.max — 18t. This defines a characteristic temperature 
Tc of the fiow to strong coupling. As it is well known, in two spatial dimensions breaking of continuous symmetries 
is impossible at T > 0. Thus the flow to strong coupling which we observe for a large part of the parameter range 
should either be interpreted as an indication of an ordered ground state of the corresponding type or as a tendency 
towards ordering at T ~ Tc for the case when an additional coupling term in the third spatial directipn, isJucluded. 

For the numerical integration of these coupled equations we apply a phase space discretizationlliJollj. The 2D 
Brillouin zone (BZ) is divided up into N elongated patches centered around N lines. Each line connects the origin 
with one of the (±7r, ±7r)-points with two straight pieces and a kink at the Umklapp surface (see Fig. ^). 

Next we discretize the coupling function. We approximate Vt(A:i, ^2, fca) by a constant for all wave vectors in the 
same patches and calculate the RG flow for the subset of interaction vertices with one wavevector representative for 
each patch. We take these wave vectors as the crossing points of the N lines with the Fermi surface (FS) . The phase 
space space integrations are performed as sums over the patches and integrations over the radial direction along 3 or 5 
lines inside each patch. Most calculations were using 48-patch systems. Calculations with 96 patches gave equivalent 
results. A typical Fermi surface with = 48 points is shown in Fig. |[ 
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The motivation for using this patch scheme is that the Fermion propagators are largest on the Fermi surface and 
decay away from it. Thus one can expect that by tracking the coupling functions for the wavevectors close to the 
Fermi surface the leading flow will be described accurately. 

Together with the flow of the interactions we calculate several static susceptibilities by coupling external fields of 
appropriate form to the electrons. During the flow these external couplings are renormalized through one-loop vertex 
corrections, as described in Rcf. In this paper, we concentrate on the (1^2 _y2- and p-wave pairing susceptibilities, 
the AF susceptibility Xs{q = (7r,7r)) and the FM susceptibility XsiQ 0). 

V. RESULTS FOR THE 2D HUBBARD MODEL AT THE VAN HOVE FILLING 

Here we describe the results for the RG flow of the coupling function for the 2D Hubbard model with initial 
interaction U = 3t and varying value for the next-nearest neighbor hopping t' . The chemical potential /i is fixedtj 
at the van Hove value /i — At' such the FS always contains the saddle points at (±7r,0) and (0, ±7r). We start the 
RG scheme at temperature Tq — 4t and integrate the coupled equations with decreasing temperature. We stop the 
flow when the largest coupling exceeds a high value larger than the bandwidth, e.g. Vr.max. — 18t, which defines the 
characteristic temperature for the flow to strong coupling, Tc- 

The overall behavior of the flow to strong coupling is shown in Fig. ^. Focusing on AF, c?-wave and FM suscepti- 
bilities, we observe three distinct parameter regions. The first, the AF regime, occurs closer to half band filling, for 
smaller absolute values of t' > —0.2t. Here the flow to strong coupling takes place at relatively high temperatures and 
the AF susceptibility is growing most strongly (see Fig. || for the Fermi surface and the flow of the susceptibilities). 
The dominant scattering processes at low T are given by AF processes between FS parts connected by wavevectors 
w (tt, tt). These can be seen in the middle plot of Fig. ^ as bright features corresponding to strong repulsive couplings 
on the line fc2 — ~ (tt, tt). 

If we increase the absolute value of t' and adjust the chemical potential such that the FS remains at the van 
Hove points, the characteristic temperature Tc for the flow to strong coupling drops continuously and for t' < ~0.2t, 
the d-wave susceptibility takes over as the leading susceptibility. In the flow of the couplings one clearly observes 
the dominant d2;2_j^2-symmetry of the pair scattering, see the diagonal features in the middle plot of Fig. ^. For 
t' « — 0.25i and a band filling slightly larger than the van Hove filling we again find a regime where the fiow of d-wave 
and AF processes is strongly coupled similar as the saddle point regime studied in Ref. |l^. 

Here we focus on the flow at the van Hove filling when we increase the absolute value of t' . The characteristic 
temperature for the fiow to strong coupling drops rapidly for t' < — 0.3t, while for t' > — 0.33i it rises again. For 
these values the fiow to strong coupling is dominated by processes with small momentum transfer (middle plot in 
Fig. 1^) and the FM susceptibility Xs{<l = 0) is by far the most divergent susceptibility at low temperatures (right 
plot in Fig. 0). The overall behavior strongly suggests a critical value for t' around t'^ « — 0.33i where the ground 
state changes from c?-wave singlet superconducting to ferromagnetic. We note that while our method allows, within 
the approximations made, to detect the instability of the Fermi liquid state against ferromagnetic fluctuations it does 
not give any information on the degree of the polarization of the ordered ground state. 

Further we stress the fundamental difference of the qualitative change in the flow between the d-wave to the FM 
regime to the crossover from the AF to the d-wave regime: in the latter case the transition is continuous and takes 
place at relatively high scales. As emphasized in Refs. |l^ and |l^, for the FS close to the saddle points and away from 
half-filling, AF and d-wave tendencies do not compete but reinforce each other on the one-loop level, therefore oue. 
finds a gradual change in the character of the flow to strong coupling when t' is varied. There are strong indicationst^ 
that the strong coupling state for a certain parameter range is not simply a symmetry-broken phase. In contrast with 
that, the transition from the dj.2_y2-wave regime to the FM regime at larger absolute values of t' is very distinct and 
- as suggested by our one-loop analysis - may be a quantum critical point of two mutually excluding tendencies. Due 
to the competition between singlet superconducting and ferromagnetic tendencies, the characteristic temperature for 
the flow to strong coupling becomes suppressed to smallest values around t'^.. The precise properties of the quantum 
critical point cannot be analyzed further for the time being because the numerical integration of the RG flow becomes 
rather time-consuming for these parameters. 

The rivalry between singlet superconducting and FM tendencies is already foreshadowed on the d-wave side at 
t' w — 0.3i. There the temperature-flow scheme yields much lower characteristic temperatures than the momentum- 
shell scheme used in Ref. |l^ which is rather insensitive against the FM tendencies arising from the van Hove points 
for the reasons mentioned in the introduction. 

The transition from singlet pairing to the ferromagnetic regime can also be seen in the flow of certain coupling 
functions which couple into both channels. For example let us take the component Vt(1, 25, 25) for a = 48 system 
(see Fig. H). The incoming wavevectors kp,i and kp,25 add up to zero, therefore the coupling corresponds to a Cooper 
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FIG. 4. Characteristic temperature Tc for the flow to strong couphng versus next nearest neighbor hopping amplitude t' from 
a 48-patch calculation (the crosses show data for 96 patches). The chemical potential is fixed at the van Hove value /i = 4t'. Tc 
is defined as temperature where the couplings reach values larger than 18t. For small \t'\, the instability is dominated by AF 
tendencies, for t' < —0.2t the d-wave tendencies prevail, and for t' < —0.33t the ferromagnetic fluctuations grow strongest. The 
criterion for the distinction between these regimes is taken as the derivative of the susceptibilities with respect to temperature 
at the scale where the couplings become larger than lOt. 
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FIG. 5. Flow to strong coupling for t' — — 0.025t and fi — —O.lt. The left plot shows the 48 points on the Fermi surface 
(points 1, 12, 13, 24 etc. are closest to the saddle points). The middle plot shows the coupling function VT(fci,fc2,fc3) at low 
r as a function of the two incoming wavevectors labeled by ki and k2 moving around the FS. fea is fixed at point 6 in the 
Brillouin zone diagonal. The colorbar denotes the value of the coupling function in units of t. The right plot shows the flow of 
the AF susceptibility (dashed line), d-wave susceptibility (solid line) and FM susceptibility (dotted line) as a function of the 
temperature T. 
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FIG. 6. Flow to strong coupling for t' = — 0.25i and /i = —t. The left plot shows the 48 points on the Fermi surface (points 1, 
12, 13, 24 etc. are closest to the saddle points). The middle plot shows the coupling function Vriki, ^2,^3) at low T as a function 
of the two incoming wavevectors labeled by fci and k2 moving around the FS. fca is fixed at point 1 close to a saddle point. The 
colorbar denotes the value of the coupling function in units of t. The right plot shows the fiow of the AF susceptibility (dashed 
line), d-wave susceptibility (solid line) and FM susceptibility (dotted line) as a function of the temperature T. 




FIG. 7. Flow to strong coupling for t' — — 0.45t and fi = —1.8t. The left plot shows the 48 points on the Fermi surface 
(points 1, 12, 13, 24 etc. are closest to the saddle points). The middle plot shows the coupling function VT{ki,k2,k3) at low 
T as a function of the two incoming wavevectors labeled by ki and k2 moving around the FS. fca is fixed at point 6 in the 
Brillouin diagonal. The colorbar denotes the value of the coupling function in units of t. The right plot shows the fiow of 
the AF susceptibility (dashed line), d-wave susceptibility (solid line) and FM susceptibility (dotted line) as a function of the 
temperature T. 
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FIG. 8. Left: Flow of the forward scattering coupling function Vt(1,25,25) for the saddle point filling and different values 
of t' . The dashed (solid) lines correspond to the AF (d-wave) regime, where Vt(1,25,25) decreases when the temperature 
T is lowered, while the dotted lines correspond to the FM regime, where Vt(1,25, 25) increases with decreasing T. Right: 
Characteristic temperature Tc for the flow to strong coupling versus t' at the van Hove filling jj, — At' . The lower curves (same 
as in Fig. ^ are obtained with the full one-loop RG, while the upper dotted line result from the RG where the particle-particle 
channel is left out. Right plot: Comparison of the flow of the FM susceptibility Xs(0) with and without particle-particle 
diagrams. For both cases, the RG equations were integrated until the couplings reach Vr.max. ~ 18f. The dashed line denotes 
the critical temperature from the generalized Stoner criterion. The vertical line at t' = — 0.32t separates the FM regime from 
the AF in the Stoner calculation and RG without particle-particle diagrams. 

process where a pair of electrons is scattered onto itself. Since wavevector kp,! is close to the saddle points it should 
rapidly diverge to —oo if the flow is towards a d-wave superconducting ground state. On the other hand Vt(1, 25, 25) 
couples into the FM channel as well, as the momentum transfer between the second incoming and the first outgoing 
wavevector is the same (in the calculation we use a very small momentum transfer). In Landau- Fermi liquid language 

Vt{1, 25, 25) corresponds to the process a°(A;F,i, — fcF,i), which should become large and positive if the flow signals a 
FM ground state. By plotting the flow of Vt(1,25,25) for the saddle point filling and different t'-values, as shown 
in Fig. ^ we can nicely observe the qualitative change in the flow when \t'\ becomes larger than the critical value 
t'^ « — 0.33i. We note that Vt(1, 25, 25) does not belong to the dominant scattering processes in the FM regime. The 
latter are given by forward scattering processes a°'{kp,k'p) with kp and k'p closer to the BZ diagonals (see Fig. |^) 
and grow more strongly towards Tc- 

One can also investigate the contributions of different one-loop diagrams to the RG flow. Here we concentrate on 
the influence of the particle-particle channel on the flow. In the left plot of Fig. ^ we compare the characteristic tem- 
peratures for the full one- loop RG with the results when the particle-particle channel is left out. For all values oft', the 
Tc is much higher, which shows clearly the importance of the particle-particle channel for a screening of the repulsive 
interaction. Without the particle-particle diagrarpa and close to t' = and t' = —Q.5t the critical temperatures from 
the RG and from the generalized Stoner-criterionEll are very similar. Again the FM susceptibility becomes dominant 
around t'^ ~ —0.33^, but without the particle-particle diagrams generating strong superconducting fluctuations, the 
characteristic temperature for the flow to strong coupling is only weakly suppressed for these parameters. 

VI. THE FLOW AWAY FROM THE VAN HOVE FILLING 

Next we investigate the flow to strong coupling for band fillings away from the van Hove filling. The flow arising for 
smaller absolute values of the next-nearest neighbor hopping il, is, ^pry similar to the flow found with the momentum- 
shell techniques which has been extensively discussed in RefsOa^E^. Therefore we now focus on the case t' < — 0.33t 
where we find a flow dominated by ferromagnetic fluctuations at the van Hove flUing. Changing the particle density 
moves the van Hove singularities away from the Fermi surface and one expects a reduction of the characteristic 
temperature for the flow to strong coupling accompanied with a decrease of the ferromagnetic tendencies. Furthermore, 
by similar arguments as for the interplay between antiferromagnetic and d-wave superconducting fluctuations when 
the FS nesting is reduced, one may presume that some kind of triplet superconducting channel might diverge at lower 
temperatures where the ferromagnetic tendencies get cut off. This expectation is confirmed by the results of the 
temperature-flow RG scheme. 

The Cooper pair scattering Vrik, —k, k') transforms according to one of the 5 irreducible representations of the 
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FIG. 9. Left plot: Flow of ferromagnetic (solid lines), px-w&ve (dashed lines) for t' — — 0.45t at fi — ~1.78t (thick lines) and 
/i = — 1.74t (thin lines). Right plot: Critical temperature Tc where the largest couplings reach 18t as function of the chemical 
potential fi. The dotted vertical lines separate the parameter region where the FM or the p-wave susceptibility grow most 
strongly. 




FIG. 10. Left panel: Fermi surface for /j, — —l.llt and t' = — 0.45t. Middle panel: snapshot of the coupling function 
V{k\, k2, fcs) with fcs = 6 fixed in the Brillouin zone diagonal. Right panel: Pair scattering V{k, —k, kg) with k varying around 
the Fermi surface at T = 1.1 • 10~^t.. For the solid line fcs is fixed at point 1 close to the saddle point (— tt, 0) while for the 
dashed line, fcs ~ ( — 1.06, —0.86) is at point 6. 



point group of the 2D square lattice. Among those there is the one-dimensional representation with d^2_j,2-symmetry 
important for smaller and the only two-dimensional representation transforming like Px and Py. In a solution 
of the BCS gap equation, components belonging to different representatfons do not mix, while e.g. within the p- 
representation components Px oc cos 6 and py oc sin 6 and the higher harmonics oc cos 30 and oc sin 30 (also called 
/-wave) will occur together. 

As the RG calculations down to low temperatures are rather time-consuming we limit ourselves to the exemplary 
case t' — — 0.45t and leave a full investigation of the parameter space to future work. As shown in the right plot of Fig. 

for t' — — 0.45i, the critical temperature drops by several orders of magnitude when we increase the particle density 
per site from the VH value (n) « 0.47 at /i = ~1.8t to (n) « 0.58 at /i = —1.7t. Further upon moving away from 
the VH filling, the growth of the FM susceptibility gets cut off and the p-wave triplet superconducting susceptibilities 
with symmetry px cx; cos 9 or py oc sin^ diverge at low temperature. Higher order harmonics cx cos 30 and sin 30 
diverge in a weaker fashion. The pair scattering at low temperatures in this parameter range is shown in Fig.|l^. One 
nicely observes that the pair scattering involving particles close to the saddle points is suppressed as the odd-parity 
nature of the p-wave pairing requires an opposite sign of the pair scattering Vrik, —k,k') e.g. for k' « (0,7r) and 
k' « (0, — tt). Comparing the temperature derivatives of ferromagnetic and superconducting susceptibilities when the 
leading couplings become larger than 12i, the transition from the ferromagnetic to the p-wave superconducting regime 
occurs at a particle density of (n) « 0.55 per site. Again, very similar to the interplay between antiferromagnetic 
and dx2-y2 superconducting fluctuations for small \t'\, on the one- loop level there is a smooth evolution of the flow to 
strong coupling from the strongly FM dominated to the predominantly p-wave dominated instability. This suggests 
the presence of a transition between the two types of ordered states as a function of the band filling provided that 
symmetry-breaking at finite temperatures becomes possible in a three-dimensional environment. From our analysis 
the transition appears to be first order, but additional interaction effects beyond our one-loop calculation could change 
the behavior. 
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angle In w.r.t. ^,7t) angle In w.r.t. (rt,7c) 

FIG. 11. Left panel: Angular variation of real (solid line) and imaginary (dashed line) part of the gap function obtained 
by the method explained in the text. Right panel: Angular variation of the gap magnitude. For these data, /i = —l.Tlt and 
t' = -0.45t. 



It is plausible to assume that the superconducting state for t' = — 0.45t will be given by a-*odeless superposition 
of the two components with symmetries Px and py, which maximizes the condensation energyl3. The direction of the 
Cooper pair spin or the (i(fc)-vector parameterizing the odd-parity gap function remains indeterminate in the absence 
of spin-orbit coupling. Therefore six different pairing symmetries d{k) oc xk^ ± yky, oc xky i ykx and oc zi^kx i 
yield the same condensation energy and additional interactions like spin orbit coupling or the spin fluctuation feedback 



mechanism are needed to select one of these representationsEj. In any case the superconducting gap function for given 
s and s' will read 

Ass'ik)^ Ao{k){kx±tky), (51) 

where the real prefactor Ao(fc) is unchanged under symmetry transformations of the lattice and takes care of the 
strong anisotropy within symmetry-related Fermi surface parts. 

Our present method does not allow a straightforward calculation of the gap function in the superconducting state. 
Nevertheless we can use the Cooper pair scattering obtained from the RG flow as effective pair potential in a mean- 
field treatment of the superconducting state. Since the Cooper scattering diverges in our approach, this will not give 
any precise information on the overall magnitude of the energy gap, but it allows us to study the symmetry and the 
angular variation of the gap function. In Fig. |l^ we show the solution of the BCS gap equation with symmetry 
given by Eq. ( |5l| ) and using the rescaled Cooper pair scattering extracted from the RG in the case /x = —l.Tlt 
and t' — —0A5t. For these data the angle is measured with respect to the (tt, tt)- point. Clearly the dips in the 



pair scattering at the saddle points mentioned above and observable in Fig. 10 lead to minima of the absolute gap 
magnitude on these FS parts. Thus the particular location of the FS close to the BZ boundary results in a highly 
anisotropic superconducting energy gap. 



VII. COMPARISON WITH OTHER APPROACHES AND CONCLUSIONS 



We have described a modified RG scheme for interacting fermion SKfterns-which uses the temperature as explicit 
scaling parameter. In comparison with the momentum-shell schemesEIEiloa the new approach has the advantage 
that it allows one to treat small-if singularities in the particle-hole channel, in our case corresponding to ferromagnetic 
tendencies, on the same basis as large-q (here antiferromagnetic) fiuctuations or superconducting correlations. 

At first sight the successive lowering of the temperature appears to be closer to the physical intuition than the 
reduction of an momentum-shell cutoff. On the other hand it should be mentioned that in general the initial conditions 
for the vertex function of the temperature-flow RG at the starting temperature do not exactly equal the vertex 
functions (to some order) of the microscopic Hamiltonian. The latter fact might complicate the comparison to other 
techniques. 

Using the temperature-flow RG scheme, we have analyzed the flow to strong coupling for the two-dimensional 
t-t' repulsive Hubbard model on the square lattice. The results obtained agree qualitatively with a various other 
calculations. 

The flow to strong coupling in the parameter region of small nex±_n£iarest neighbor hopping t' shows strong similar- 
ities to the RG flow known from the momentum-shell techniquealilo'Ej. Again AF and d-wave tendencies are coupled 
to a large extent and the FS curvature determines which of the two types of fluctuations is strongest: when we move 
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away from perfect nesting, the AF processes are cut off at low temperatures and the dx2_y2--wave Cooper scattering 
processes can diverge. 

For larger absolute values of t' beyond a critical t'^ « — 0.33i we find that the flow to strong coupling is dorniiiated. 
by ferromagnetic fluctuations. This regime has not been found with the momentum-shell schemes used recentlyEJo'E£l 
for the reasons explained in the introduction. _ 

There are several other approaches which predict ferromagnetic ground states for larger \t'\. Lin and HirschE^ 
analyzed the Hartree-Fock phase diagram for the t-t' Hubbard model and found the ferromagnetic state to be stable 
against the antiferromagnetic state for t' < — 0.324i. They do not find a d-wave phase for t' > — 0.324i because the 
bare repulsive Hubbard interaction does not contain any attractive component in the pair scattering channel. The 
critical t'-value from Hartree-Fock is practically equivalent to the critical t'^ ~ — 0.33t from the RG treatment. On 
the other hand in Hartree-Fock one would find a first order transition at t'^, while in our treatment the characteristic 
temperature gets suppressed to zero at the quantum critical point. Irkhin et al.E3 applied a parquet scheme using the 
same set of one-loop diagrams as the RG described above. With a simplified dispersion they found indications for 
ferromagnetism close to t' — —0.5t [and. emphasized the strong deviations of the critical temperatures from the simple 
Stoner calculation. Hlubina et al.EfrES analyzed the question of ferromagnetism by applying quantum Montecarlo 
and T-matrix approximation (TMA). The latter method gave a ferromagnetic state beyond a critical t' « — 0.43i. 
The Kanamori-screening in the particle-particle channel of the TMA is included in our temperature-flow scheme, 
which gives a larger window for a ferromagnetic ground state. This suggests that the straightforward TMA might 
overestimate the screening effects in the given example. We repeat however that defining the initial condition for two- 
point and four-point vertex functions in the temppature-flow RG is not fully equivalent to using a specific microscopic 
Hamiltonian in the perturbative scheme. Aracheao found spin-polarized ground states in exact diagonalization studies 
of a 4 X 4 t-t' cluster for t' = -a|it4»d t' = -0.6t at densities 0.5 and 0.375. 

For several classes of modelstS^cS-ES large-spin ground states can be proved exactly. Among those, the so-called 
fiat band ferromagnetism appears to have the closest relation to the case of the t-t' Hubbard model. The typical 
feature of these models is a large density of states at the bottom of the band which gets pushed further below the 
Fermi energy by polarizing the system. In the 2D t-t' Hubbard model for t' — 0.5t the density of states is peaked 
at the bottom of the band. Nonetheless from the perspective of our RG calculation the large density of states at the 
FS appears to be essential for the FM tendencies in the t-t' Hubbard model, as the latter are strongly reduced when 
one raises the Fermi level above the van Hove energy. Moreover, if we cut out the bottom of the band by hand and 
discard all modes with eg < — O.lt, the ferromagnetic tendencies do increase rather than decrease, presumably due to 
the reduced screening of the onsite repulsion. Thus the peak at the bottom of the band is not essential for our finding 
of ferromagnetic tendencies. The sensitivity of the ferromagnetic regime to the location of the FS points to another 
effect which is not included in the present calculation. It is probable that the inclusion of sclfenergy effects in the 
fiow will shift the Fermi surface and alter the low energy density of states such that the actual ferromagnetic regime 
might be changed. A calculation of the selfenergy at the FS from the fiow of the interactions shows however that the 
real part of the selfenergy becomes comparable to the temperature only when the couplings have grown larger than 
the bandwidth to Vr.max. « 12 — 13t. 

In analogy with the interplay between AF and d-wave fluctuations close to perfect nesting at smaller t', there exists 
a parameter region for larger \t'\ where the FM tendencies get cut off and the flow to strong coupling is dominated 
by Px- and p^-wave superconducting correlations. The energy gap suggested by the RG flow of the Cooper pair 
scattering is highly anisotropic and exhibits paipima. on the FS parts close to the saddle points. Similar scenarios 
have already been proposed by several authoraHJ'OEll, our approach allows a more systematic analysis of the crossover 
from the FM to p-wave regime and the possible order parameter symmetries. We note that the superconducting gap 
functions suggested by the RG flow (x Px -\- ipy are among the candidates under debate for the superconducting state 
of Sr2Ru04c3. In this quasi-two-dimensional system one of the three Fermi surfaces has a similar shape as the cases 
studied above. It should be interesting to investigate the RG flow and the superconducting properties suggested by 
the resulting pair scattering in a more realistic three-band model. 

Another future application of the modified RG scheme could be to investigate the interplay of screening of long-range 
forces with other tendencies. Further, other potential instabilities in forward-scattering channel could be analyzed, 
e.g. Labbe-Etiedel or Pomeranchuk instabilities with spontaneous deformations of the FS, as proposed by Halboth 
and MetzneitJ. 
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